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Abstract 

We study the effect of thermal fluctuations of the s-wave order parameter 
of a quasi two dimensional superconductor on the nuclear spin relaxation 
rate near the transition temperature Tq- We consider both the effects of the 
amplitude fluctuations and the Berezinskii-Kosterlitz-Thouless (BKT) phase 
fluctuations in weakly coupled layered superconductors. In the treatment of 
the amplitude fluctuations we employ the Gaussian approximation and eval- 
uate the longitudinal relaxation rate T^ 1 for a clean s-wave superconductor, 
with and without pair breaking effects, using the static pair fluctuation prop- 
agator T>. The increase in T-j" due to pair breaking in T> is overcompensated 
by the decrease arising from the single particle Green's functions. The re- 
sult is a strong effect on T^ 1 for even a small amount of pair breaking. The 
phase fluctuations are described in terms of dynamical BKT excitations in 
the form of pancake vortex-antivortex (VA) pairs. We calculate the effect 
of the magnetic field fluctuations caused by the translational motion of VA 
excitations on T{~ and on the transverse relaxation rate TrT 1 on both sides 
of the BKT transition temperature Tbkt < Tq. The results for the NQR 
relaxation rates depend strongly on the diffusion constant D that governs the 
motion of free and bound vortices as well as the annihilation of VA pairs. We 



discuss the relaxation rates for real multilayer systems where D can be small 
and thus increase the lifetime of a VA pair, leading to an enhancement of the 
rates. We also discuss in some detail the experimental feasibility of observing 
the effects of amplitude fluctuations in layered s-wave superconductors such 
a the dichalcogenides and the effects of phase fluctuations in s- or d-wave 
superconductors such as the layered cuprates. 



I. INTRODUCTION 

The most common NMR and NQR experiments on high-Tc and other quasi two di- 
mensional superconductors concern the Knight shift and the longitudinal and transverse 
relaxation rates. Both of these experiments explore the low-frequency spin dynamics of elec- 
trons and holes in normal metals and superconductors. The relaxation rates are caused by 
the time dependence of the fluctuating magnetic fields. In superconductors these field fluc- 
tuations originate from electronic quasiparticle excitations and from the motion of magnetic 
vortices. In the high-Tc cuprates the quadrupolar Cu spin-lattice relaxation can be due to 
the transitions between quadrupolar states of the Cu nuclei caused by the interaction of the 
nuclear spins with the time dependent magnetic fields of the vortices. The 2D cuprates such 
as Bi-2212 consist of Cu02 layers with a very small interlayer hopping matrix elememt t^.Q 
The magnetic field fluctuations near the real superconducting transition temperature Tc , 
where the resistivity goes to zero and the long-range order is established by Josephson phase 
coupling, are caused by both the quasiparticle excitations of the normal and superconduct- 
ing states and by the spontaneous excitation of thermal vortex-antivortex (VA) pairs. These 
vortex excitations occur close to the Berezinskii-Kosterlitz-Thouless (BKT) transition tem- 
perature of a 2D layer, Xbkt < Tc , cf. Fig. Whereas the long- wavelength Gaussian 
fluctuations of the quasiparticle excitations consist of amplitude and long-wavelength phase 
fluctuations of the complex order parameter, the VA fluctuations are primarily phase fluc- 
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tuations of the order parameter. The BKT scenario also applies to 2d layers that are weakly 
coupled by electromagnetic and Josephson coupling effects. Around the real transition 
temperature Tq there exists the narrow Ginzburg regime of critical fluctuations where the 
Josephson and Coulomb couplings into the third dimension begin to establish the phase 
coherence between neighboring layers. This is representated by the shaded area around Tq 
in Fig. 1. According to the detailed calculations of the specific heat fluctuations of Ramallo 
and Vidali and their comparison with experimental results, the three-dimensional critical 
behavior is restricted to a rather narrow range in the cuprates, t = |T — Tc\ < 10~ 2 , so 
that a wide region of Gaussian fluctuations of the order of 10K exists above To- The same 
considerations should be be valid for the dichalcogenides which are layered systems exhibit- 
ing s-wave pairing. The BKT vortex-antivortex fluctuations exist above Tq below the mean 
field transition temperature T C o (cf. Fig. 1). At T cr > T c , the inter layer phase coherence 
becomes so weak that the Josephson vortices proliferate and the 3D phase coupling ceases 
to exist. 

In Sec. II we study the effect of Gaussian fluctuations on the spin-lattice relaxation rate 
of s-wave superconductors near the transition temperature for layered systems. The NMR 
spin relaxation rate T^" 1 has been investigated intensively in the high-Tc superconductors. 
At present however the general opinion seems to be that only the electron-doped high-Tc 
superconductors may exhibit s-wave pairing.i'0 Other layered compounds with s-wave pair- 
ing are the layered transition metal dichalcogenides. Long wavelength thermal fluctuations 
of the superconducting order parameter, i.e., Gaussian fluctuations, are expected to play a 
pertinent role in layered compounds for several reasons; in particular, the quasi 2D structure 
of these metallic systems, the small coherence length £o m the layers (cuprates), and the 
large penetration depth A for fields parallel to the layers. The effect of such fluctuations 
on T{ 1 in these and similar systems has been the subject of numerous papers. Most of 
the previous authorsSSEHJIl] included non-magnetic impurity scattering and also set the 
small external NMR frequency u equal to zero. In such theories it is not clear whether the 
correct clean limit is obtained. For example, in a clean system without pair breaking, we 
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find 1/Ti oc 1/u , for T — > Tc . To elucidate this limit we consider a clean system from the 
outset. This may be a reasonable first approximation for systems like the cuprates where £o 
is large compared with the mean free path £. 

Certainly for the cuprates, and many other layered systems as well, it is important to 
include the pair breaking effect of inelastic scattering due to exchange of low energy Bosons, 
spin fluctuations or phonons, for example. To do this correctly would require solution 
of the strong-coupling Eliashberg equations taking full account of the energy dependence 
of the pairing interaction and the single-particle self energy. To date no such theory has 
appeared in the literature with respect to NMR in layered systems for s-wave pairing.il In 
the previous weak-coupling theories inelastic scattering has usually been accounted for by 
introducing a cutoff in the pair fluctuation propagator or a constant lifetime in the single 
particle Green's function. We also do a weak-coupling calculation and simulate the effect 
of inelastic scattering through a pair breaking parameter in analogy to magnetic impurity 
scattering. The large qualitative and quantitative effect of pair breaking indicates that 
the results of a weak coupling calculation should be viewed with caution until they can be 
confirmed with a strong coupling Eliashberg calculation. Since the inclusion of pair breaking 
in the clean case qualitatively changes the results we give some of the details both with and 
without pair breaking. 

In Sec. Ill we study the effect of topological intra-layer phase fluctuations (vortices) on 
the NQR spin-lattice and spin-spin relaxation rates T{~ and T 2 _1 . In contrast to the Gaus- 
sian fluctuations of Sec. II, which contain the long-wavelength fluctuations of the amplitude 
of the order parameter, the topological fluctuations are, primarily, the 2D singular phase 
fluctuations of the superconducting order parameter. Experimentally, the relaxation rates 
in the cuprates can be obtained in zero magnetic field from NQR spin-echo experiments 
on 63 Cu and 65 Cu nuclei. We consider NQR to allow direct application of the zero-field 
BKT theory. The translational motion of vortices and antivortices and the corresponding 
time dependent magnetic fields not only affect the nuclear spin relaxation but also the NQR 
line width and for this reason we calculate, in addition to T-f 1 , the spin-phase relaxation 
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rate T 2 ~ . The distribution of VA pairs is determined by the singular parts of the otherwise 
smooth phase fields of the complex order parameter. The topological VA excitations reduce 
the transition temperature of a single layer to the Berezinskii-Kosterlitz-Thouless tempera- 
ture Tbkt j above which the bound VA pairs begin to break up into free pancake vortices 
and antivortices. We study, in the temperature regime above Tbkt, the magnetic field cor- 
relation functions that determine T± and T 2 _1 . The magnetic field fluctuations are caused 
by the translational motion of the thermally excited vortices and antivortices. The diffusion 
constant of the vortices and antivortices, D(T), plays a crucial role for the magnitudes of 
the relaxation rates; it determines both the free vortex and antivortex motions and the re- 
combination time of a VA pair. We find that D must be sufficiently small in order to get 
vortex relaxation rates Tf"^ that are comparable to those caused by quasiparticle relaxation. 
The latter cause relaxation rates of the order of 10 2 to 10 3 s -1 near Tc of the cuprates. The 
effect of the time- varying magnetic fields due to the VA fluctuations on T 2 _1 is also of some 
interest for the following reason: Of the two mechanisms causing the transverse relaxation 
in the cuprates, namely the indirect nuclear spin-spin coupling and the spin-lattice relax- 
ation, we expect that only the second contribution is affected by the VA fluctuations. These 
fluctuations lead to an exponential decay of the time dependent transverse magnetization 
of the 63 Cu or 65 Cu nuclear spins.S On the other hand, the individual 63 Cu - 63 Cu indirect 
coupling leads to a Gaussian time decay and is of such short range, r <C £ ,0 that this 
coupling remains altogether unaffected by the transition into the superconducting state. 

In Sec. IV we discuss the experimental situation in some detail with respect to the 
possibility of observing the effect of superconducting fluctuations on the longitudinal and 
transverse relaxation rates of the nuclear spins in layered systems such as the cuprates 
with s-wave pairing and conventional layered superconductors such as the transition metal 
dichalcogenides. 
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II. EFFECT OF PAIRING FLUCTUATIONS ON Tr 1 



A. General Formalism 



The NMR spin relaxation rate is given by 



1 _ EAq W(q 1 oO ) (1) 



where x is the transverse (spin-flip) susceptibility, u is the external frequency, and A q is 
determined by the fine structure constants and should be large for the q's that couple to large 
values of Imx(q, u). Here we assume that A q is a constant. The susceptibility diagrams 
through first order in the Cooper pair fluctuation propagator, V, are shown in Fig. 2. 
The zeroth order contribution, x yields Korringa's law in the normal phase and, in the 
absence of pair-breaking, the Hebel-Slichter peak below Tc for s-wave pairing. The leading 
fluctuation contributions are given by the diagrams xff and Xgg ■ We follow the notation 
of Maniv and Alexander^ , who considered the clean case in 3D, and write 

2 Imx(q, u) = ImxH, (2) 
q 

where 

ImxM = ImxFF^) + ImxGcO) • (3) 

The first term represents the Maki-Thompson (MT) diagram,!!!! Xff = Xmt, which also 
has a long history in the calculation of the fluctuation conductivity. The second term 
accounts for the self energy, or density of states (DOS), effect that is caused by the renor- 
malization of the normal state Green's functions by superconducting fluctuations. This 
contribution is often denoted by Xdos — Xgg The second-order (Aslamazov-Larkin) 

diagram does not contribute to transverse susceptibility. 

We define the static fluctuation propagator T> in terms of the particle-particle t-matrix 



as 



X>(k) = -fc B Tt(k,0), (4) 



where 

d 3 k' 

£ _1 (k, iu m ) = V' 1 - k B T J2 TTls^^' iuJ n)G{<\ - k', iv m - iu n ) 7 . (5) 

Here u; n = u n (T) = (2n + l)Trk B T, k is the total momentum of the Cooper pair, 7 is the 
impurity vertex, and we have taken the pair interaction V to be constant. P(/c)diverges 
when T approaches Tq from above for k equal to zero, signaling the transition to the su- 
perconducting phase and defining Tq- In order to describe static Gaussian fluctuations, 
the t-matrix is calculated for small total momentum with zero order Green's functions. A 
standard weak coupling calculation!! yields for a clean system (7 = 1) 

T)-i(v\ = f^L] J V kBTc - V kBT i- (fe^l! v keT 1 (6) 

™* T ' uk) \ u »F°)\ Jfo m^i 12 kcoi 8 J • 1 } 

where = ma 2 /(2irTi 2 ) is the 2D density of states per spin and Vp is the Fermi velocity. 
All the sums are cut off at the constant BCS cutoff uc- The cutoff on n, n c = — |, 
depends on T and is responsible for the important difference in the first two sums in Eq. 
(|D .H Performing the sums leads to 

_ k B T/N 

U{k) ~ \n(T/T c ) + C 2 k 2 ' [ V 

where £ is related to the BCS coherence length £ = £(T = 0), 



/ 7C(3) 7w F f ~ /T1 , fi v 
* = V 48 nk B T=^ Tc/T ^ (8) 

In the early conductivity calculations it was first observed that the Maki diagram has anoma- 
lous properties: In the presence of non-magnetic impurity scattering, the impurity vertex 
corrections at each end of the fluctuation propagator T> in the polarization bubble lead to 
a divergence at all temperatures in 2D. This divergence could be removed by including a 
pair breaking parameter.00 We emphasize that we work in the clean limit where there is 
no impurity vertex in the first place and thus this "Maki-divergence" does not occur. As we 
will see however, this contribution is also rather pathological in a clean system in 2D: We 
find that it diverges at all temperatures when u goes to zero. This divergence is removed 
(except right at T = Tc) by inclusion of pair breaking. 



The Maki-Thompson contribution to Tf , xff, was first evaluated in 2D by Kuboki 
and Fukuyamalll in the dirty limit for u = 0. They found a logarithmic divergence on 
approaching Tq from above with the form \n(T/t) where T is a pair breaking parameter. 
This presented a problem in the high-T^ superconductors because no corresponding peak 
has been observed there. An attempt to remedy this situation was made by Heym0 who 
generalized Kuboki and Fukuyamas work by including the frequency dependence of P. He 
found that dynamic fluctuations becomes important for T/Tq > 1.05 and can lead to a 
significant correction to the T-dependence of 1/Ti for these temperatures, in particular for 
strong pair breaking. The large peak for T — > Tc remains however. We will see later that 
pair breaking can strongly reduce this peak. 

We follow the notation of Maniv and Alexander^! and write in 2D 

XFF(M u ) = J2c(k,t{i u )v(k), (9) 

k 

xgg^u) = J2 A (^^Mk) , (io) 

k 

where the contributions of the polarization diagrams are given by 

C(k,iQ u ) = -k B TY^ G(pi,iu n )G(k - pi, -iu n ) 

n Pi,P2 

x G(p 2 , ioJ n + zO^)G(k - p 2 , -ioJ n - iVt v ) (11) 

and 

A(k,iQ u ) = 2k B TY J C 2 (pi, iu n )G(k - pi, -iu n )G(p 2 , iv n + ift v ) ■ (12) 

" Pi,P2 

The G's are taken to be zero order propagators of the form (for the case of no pair breaking) 
G(p, iuj n ) = {ioj n — £ p ) _1 where we assume free particles in 2D, £ p = p 2 /2m — Ep ■ Since 
T>(k) is peaked at small k, we approxoimate £k- P ~ (, P — vpkcoscj) where = Z(k, p). The 
factor 2 in ([T2| ) accounts for the two diagrams of the GG type. 



S 



B. Clean System without Pair Breaking 



We consider first the Maki-Thompson diagram. Performing the integrals over £ pi and 
£ P2 in Eq. (|TT|) we have 



C(k, ifl u ) = Nq7i 2 T'^2 sgn{uj n )sgn(uj n + Vt u ) 



x 



lo 2n Jo 2ir (iu n — v ik/2)(iu n + iVt u — v 2 k/2) 

Due to the sgn functions it is mathematically convenient to follow Maniv and Alexanderlll 
and split the contribution into two parts: 



(2) 



(14) 



where 



27T JO 27T 



n j -(£i,£ 2 ,^) 



(15) 



with £j = kvF cos0j/2, j = 1,2 and 

+oo 



(iu n - si)(iu n - e 2 + i£l v ) 



U 2 (e 1 , e 2 , iQ v ) = -2n 2 k B T 



E" 1 

Z^n=0 



1 



1 



(iwn~ei)(iu) n -£2+i£lv) 



v > 



f < 



(16) 



(17) 



For a clean system this decomposition is done merely for mathematical convenience When 
impurity scattering is present however, the singular q dependence due to the impurity ver- 
tices at the ends of the T> propagator (the "Maki divergence" , again) makes it necessary to 
split xff i n a different manner into a "regular" and an "anomalous" part.0 Performing the 
sum over n in Eq. fll6| ) and analytically continuing iQ u to the real external frequency uo , 
we find, to leading order in to 



T T F ^2 



x , f l ~ m/x dy (-l)f'(e c xy/T) 
V(x)dx j = ., ^;= J = , (18) 



i-(f + v)' 



V(x) 



V(k) 



(19) 



k B T v "' x 2 + t(T/C) 2 ' 

where we have defined x = k/kF, ^ = hu/2E F , e c = EpjksTc , T = T/T c , C — ^f£o an d 
t = InT ~ (T — Tc)/Tc- Here /' is the derivative of the Fermi function and A is a cutoff 
which we take as 0(1). A larger cutoff does not change the results significantly. Note that, 
since — /' is positive, Imx^p is a positive contribution, proportional to Nq. 

After analytic continuation and expansion to first order in u we find from Eqs. ([T^) and 

(0) 



/ a \ 2 %uj N 








\To) E~ f ~2 




InT. 


1 



X 



277 



2tt 



oo 



xT>(x)dx 



2 2 5 
Pn^\^2 ~ Pn 



(20) 



where fa = e c x cos <pi/T , p n = (2n + 1)tt, and several terms that vanish after the angular 

(2) 

integrations have been omitted. Imxpp is negative. 

We turn now to the calculation of the DOS contribution xgg an d first evaluate the sum 
over n in Eq. ([12]) as a contour integration and then carry out the analytic continuation 
ifl nu — > uj + i8 . Of the several resulting terms only the following yields, after the angular 
integrations, a non-vanishing contribution to first order in u : 

-/(6R 2 



/+oo 
'OO 



(21) 



Pl - -oo (& - £l - o> - 2<$) 2 (6 +£i-v-k-u-id) ' 
The ^2 integral can be expressed as a sum of the residues at the poles of the Fermi function: 

1 



A{k,u) = 8<K 2 tN k B Tj2 E 



(22) 



Pl n=-1 {iuj n - & - w) 2 (i^ n + - v • k - u) ' 

Expanding to first order in u and performing the £i integration we obtain the negative 
contribution 



Ihixgg^) 



To E~ F ^2 



TFT 



A «r/2 

xV{x)dx \ i 
o Jo 



•E 

n=0 



p 2 - 3(xe c cos0/T) s 



p 2 + (xe c cos0/T) 2 



(23) 
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In Fig. 3 we plot the three contributions to Im^ as functions of the reduced Temperature 
T = T/Tc . The positive contribution Imx^pp is seen to dominate strongly. 

It is interesting to consider these results in the limit of T — > Tq from above. The 
experimental external frequency to is very small but we leave it finite since ImxFF depends 
on the ratio t/u. In Eqs. (p0|) and ( |23|) the function xV(x) oc x/(x 2 + tT /( 2 ) is peaked at 
small x for t small. The sum over n is only weakly dependent on x for small x and it is a 
good first approximation to set x = within the sum. The t dependence then arises solely 
from the ^-integration over the fluctuation propagator and one easily finds 

l^ Klnt , !^ Kln( , fori ^ . ,24, 

UJ UJ 



These limiting functions are independent of u . In Eq. (|18| ) for \mxpp the limiting behavior 
is not just due to the fluctuation propagator but is strongly affected by the integral over y. 
We find the following result: 

Ihiy (1) f -77 In (^-) for ^- -> 
- for ^ ^ 



UJ 



For comparison we give the corresponding results for small t and uj for the clean 3D case:l 
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oc In — , oc c 2 , oc -c 3 , (26) 



UJ \t + C\UJ l J UJ " UJ 

where Cj, C2, and C3 are positive constants. In comparison to 3D, the results in 2D are 
rather pathological. In 3D the only divergence occurs in lm.Xp p/uJ and then only when both 
t and uj go to zero. In 2D, Imx^/cu and Imxcc/ 1 ^ diverge only for t — > , while Irnx^/u; 
diverges as uj — ► for all t ! This strange behavior is presumably unphysical and is removed 
by a small amount of pair breaking. As seen from Eq. (p5[) , the limiting value of Imxpl?/ 00 
depends on the ratio uj/t . As t — > there is a "crossover" near t = uj/Ep which is probably 
too close to t = to be experimentally observable. In any case, for finite uj there is no 
divergence in Imxpp/uJ for t — > . 

It is also interesting to compare our Eq. ( p5| ) with the corresponding result of Randeria 
and Varlamov (RV)0 who give a result for the "ultra-clean" limit of a theory that includes 
elastic scattering from the beginning. Their Eq.(16) reads 

11 



ImXFF oc — y= In (t c T\/tj for T c ry/t 3> 1 and uj = , (27) 

where r is the lifetime for elastic scattering. Since this equation does not contain their pair 
breaking parameter S, it should be valid for the case of no pair breaking. Equations ( p5|) and 
(p7|) are similar in several respects. In both cases there is a prefactor 1/y/t multiplied by the 
In of a large number proportional to \fi . Since RV set uj = before the calculation while 
we have 1/r = , an exact comparison of the results is not possible. In the exact clean limit 
Eq. (p5]) seems preferable because uj , although very small, is a well defined experimental 
quantity while r in Eq. fl2~7|) is not well known and , more importantly, r should not even 
appear in the exact result for a clean system. Also, the limit t — > is only possible if uj 
remains non-zero. 

C. Clean system with pair breaking 

With "pair breaking" we mean essentially the effects of inelastic scattering which, in the 
high Tq superconductors, for example, could be due to spin fluctuations and phonons. In 
order to obtain a rough estimate of this effect we want to simulate it in a simple manner 
within a weak coupling theory. Our procedure is equivalent to adding a constant inelastic 
scattering rate to the single particle propagator as done by other authors.0 We attempt to 
justify this physically by assuming that the effect on the fluctuation propagator T> is simi- 
lar to the well known pair breaking effect of scattering by magnetic impurities.^ Impurity 
scattering affects T> in two ways: self energy corrections to the single particle propagators of 
the t-matrix ladder and vertex corrections to the pair interaction. For non-magnetic impu- 
rities these two contributions cancel for s-wave pairing (Anderson's theorem). For magnetic 
impurities there is no cancellation and the transition temperature is strongly reduced. For 
our purposes it is sufficient to retain only the self energy corrections and to assume that the 
vertex correction is included in the effective pairing interaction. The Green's functions in 
Eq. (||) thus have the form 
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G(k,iu; n ) = - (28) 

tw n — ^ k 

with 

u n = u n + Tsgn(u; n ) , (29) 

where T = Ti/Zr^ is the phenomeno logical pair breaking parameter destroying the phase 
coherence between the Cooper pairs. Eq. (§) for T> is now replaced by 



_ ( N o 7[ \ J \p k B T C o _ k B T _ (hv F k) 2 ^ k B T \ 



where Tc is the transition temperature in the absence of pair breaking. Although we report 
here only results for V constant, in general V will be a function of u n . For comparison, for 
T constant the sums can be carried out analytically and Eq. (p0|) can be expressed in terms 
of the digamma and tetragamma functions in the notation of Eq. (|T9|), as 



Tcq/ V2v \2 



Vl4C(3)T 2 , 



* (2) Q + «), (31) 



where 



« = — • ( 32 ) 

The transition temperature in the presence of pair breaking, Tq , is obtained by setting 
ViO)- 1 equal to zero yielding the Abrikosov-Gorkov equation£3 Since we assume the pair 
breaking is due to inelastic scattering, Tco is not experimentally accessible. Thinking for the 
moment of the high Tc superconductors, in our numerical calculations we take Tc = 100 K 
as experimentally given and Tco wm De determined by the damping parameter T . Although 
Tc < Tco j as expected, T> as a function of k can be modified by the pair breaking in such 
a way that, neglecting the effect of pair breaking in the rest of the susceptibility diagrams, 
the fluctuation contributions to T{ 1 are increased. To elucidate this behavoir we show T> 
as a function of x — k/kp in Fig. 4 for T = 1.03 Tq and several values of «c — aT. 
The figure shows that the maximum and width of T>(k) increase with T. The fact that the 
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fluctuation propagator T> increases with increasing pair breaking may seem counter-intuitive 
if one assumes, in analogy to the case of magnetic impurities where Tco is known, that: pair 
breaking — > reduced tendency to superconductivity — ► weaker fluctuations. In that case, 
for T fixed relative to Tco, the fluctuations weaken with increasing T because Tc, and thus 
the divergence of T>, are moving away from the reference point T. In our case, however, T is 
fixed relative to Tc and Tco moves away with increasing T. Thus the fluctuations at T would 
be expected to be rather independent of T . That T> actually increases arises mathematically 
from the dependence of the In term in Eq. (|3lD on T in this case. The effect of pairbreaking 
in the single particle propagators however leads to a decrease of the susceptibility diagrams 
which usually dominates over the effect arising from T>. 

We consider first x^ff si nce it is the dominant contribution without pair breaking and its 
limiting t-dependence is qualitatively changed by the addition of pair breaking. Equations 
(|i~5|) and (|16|) are still valid if the frequencies are renormalized according to Eq. (p9| ) . The 
frequency sum and then the angular integrations can be carried out exactly with the result: 



a 



Rf~1T 



(33) 



e fA- r+oo 

^= .„ t=tz / dxxV(x) duf(27ru)Y(u,r],a)[Y(u + tt,T],a)-Y(u-n,r),a)], 

lD7T 4 i \ l JO J-oo 



(34) 



where f is the Fermi function, r\ = e c x/2iiT, Q = Tiuj/2-KkBT, and Y(a,b,c) = 
tt^2 (b 2 - a 2 + c 2 + v^) IN 1/2 , with N = (b 2 - a 2 + c 2 ) 2 + 4a 2 c 2 . We have computed 
Jr numerically and the result as a function of the pair breaking parameter ac is shown in 
Fig. 5 for several values of the external frequency, Q = %uj/2Ef ■ Note the very large and 
rapid change at small pair breaking for the experimentally relevant frequency 
= 5x 1(T 8 . 

(2) 

Caculation of Imx FF (u>) for finite T yields 



a\ Two No 

To fTf^ 



ttT 



A r'l- 

xT>(x)dx 



2tt 
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X 



2vra c (/if + pl + J + (/i 2 - fii)fiip n+ 



to Oh ~ ^i) 2 + (^ac) 2 \ (p 2 + + /i?) 2 

[27Ta C (// 2 + p£_) + (A*2 - V\)^2Vn- 



(Pl- +^ 



(35) 



where p„± = p n ± 2-nac/T . A similar calculation leads to the result that Imx GG) r(^) is 
simply given by Eq. fl2"3|) with the replacement p n — >■ p n + 2irac/T . 

The limiting behavior of Inr\; GGi r(^) and Imx^ r (u;) for t = (T — T c )/T c — > is still 
given by Eq. (^). Similarly, for finite T , the T-dependence of Imx^ r (u;) now also arises 
primarily from the integral of xD{x) in Eq. (|3~4D and one finds Imx^ r (a;) oc — lnt, for 
t — > 0. Thus, in the presence of pair breaking, the magnitudes of all three contributions 
have the same limiting Temperature dependence. 

In Fig. 6 we plot the fluctuation contribution to [l/(TjT)]pT , Imx r (u;), given by equa- 
tions (|lj) and (fj), vs etc- Note that ImxpFr increases in magnitude with increasing a c until 
it dominates over Imx^p leading to a change of sign of the total contribution (dotted line) 
for the pair breaking parameter otc near 0.05. This effect also occurs in the presence of weak 
elastic scattering.0 In Fig. 7a we plot the total contribution vsT/Tc for a range of otc from 

to 0.1. The detailed plots in Figs. 7b and 7c, for small and large ac, show the dominance 

(2) 

of Imxppr f° r large pair breaking yielding a negative divergence. We point out that, in the 
presense of inelastic scattering (pair breaking), our results here in the exact clean limit are 
qualitatively similar to those with weak elastic scattering.0 

To summarize briefly, we have shown that the Maki diagram in 2D is pathological in 
the exact clean limit but the divergences occur for a different reason than in the usual dirty 
limit. In the presence of pair breaking due to inelastic scattering, which of course will always 
be present to some extent, reasonable results are obtained. A small amount of pair breaking 
also strongly reduces the increase in T{ 1 as Tq is approached. We have also seen that pair 
breaking in the fluctuation propagator V can affect Xf 1 quite differently than pair breaking 
in the Green's functions in the remainder of the diagram. In ImxpFri 00 ) 5 f° r example, the 
increase in due to pair breaking in T> is more than compensated by the decrease arising 

15 



from the Green's functions. 

The strong effect of even a very small amount of pair breaking within our simple weak 
coupling model underlines the need for a strong coupling Eliashberg calculation including 
inelastic scattering due to Boson exchange before quantitative comparison with experiment 
can be attempted for specific systems. 

III. EFFECT OF BKT VORTEX- ANTIVORTEX FLUCTUATIONS ON Tf 1 

AND TjT 1 

Up to now we have discussed the effect of Gaussian fluctuations of the order parameter on 
the NMR relaxation rate T{~ near the mean field transition temperature Tcq of a 2D super- 
conductor. The weakly coupled 2D layers of high-Tc superconductors promote the formation 
of topological excitations in the form of pancake vortex-antivortex (VA) pairs associated with 
the singular part of the phase field of the complex order parameter. Neglecting Josephson 
coupling, these VA excitations reduce the transition temperature to a Berezinskii-Kosterlitz- 
ThoulessS (BKT) transition temperature Tbkt, above which the bound pairs start to break 
up into free vortices and antivortices. We now proceed to study the effect of magnetic-field 
fluctuations caused by the translational motions of vortices and antivortices on Tf 1 and T 2 _1 
in the vicinity of Tbkt- 



A. Vortex Fluctuations in BKT Theory 

In thin superconducting films, BKT theory0 predicts that below a transition temper- 
ature Tbkt spontaneously created pancake vortices and antivortices (VA's) are bound in 
pairs with zero total magnetic flux and do not destroy the off-diagonal quasi-long-range 
order. Above Tbkt ; the large pairs break up into free vortices and antivortices, which 
are responsible for the dissipation of electrical currents, and quasi-long-range order is lost. 
Around Tbkt there is a vortex-antivortex fluctuation regime. The time and distance be- 
havior of these fluctuations affect the transport properties, e.g., BKT behavior is clearly 
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seen at microwave frequencies! and in the DC current-voltage characteristics. This pict 
is essentially unchanged in layered superconductors if the Josephson coupling is ignored.!!'! 
(Under "Josephson coupling" one understands all interlayer pair transitions contributing to 
the Josephson current.) 

In weakly coupled high-T c superconductors such as E^S^CaC^Og+tj (Bi-2212) the re- 
sults are changed, since the pancake vortices in the layers become connected by Josephson 
vortices between the layersS'i'ii The Josephson vortices lead to a linear term in the interac- 
tion of pancakes connected in this way. The general picture is the following!'!: The VA pairs 
start to unbind at a BKT temperature Tbkt, but the linear interaction leads to confinement 
of the pairs at a length scale A, called the Josephson length. The BKT renormalization is 
cut off at pair sizes of the order of A. Very close to Tbkt, where the BKT correlation length 
£bkt exceeds A, the interlayer coupling becomes important and the system shows three- 
dimensional critical behavior.il The true transition takes place at a temperature Tq- At a 
higher temperature, where the BKT correlation length £bkt(T) falls below A, 2D fluctua- 
tions become important again and remain essential up to Tco, where the local condensation 
energy vanishes. Critical fluctuations shift Tqo slighty downwards from the mean field value. 
Around Tc there is a narrow 3D critical regime in which fluctuating vortex loops through 
more than a single layer are crucial. We are considering weakly coupled high-Tc supercon- 
ductors in the sense that this 3D critical region is much narrower than the 2D fluctuation 
regime, which is well described by BKT theory. We therefore neglect the Josephson coupling 
and restrict ourselves to the 2D fluctuation region. 

In the following we are concerned with the effect of BKT VA fluctuations on the spin- 
lattice relaxation rate Tf 1 and the spin-spin relaxation rate T 2 _1 . Up to now, the experiments 
on BKT fluctuations consist mainly of flux-noise measurement^^ with frequencies of the 
order of uj ~ 10 4 s _1 , whereas NMR and NQR spectroscopy are characterized by frequencies 
of order 10 7 to 10 8 s _1 . 
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B. Vortex Relaxation Mechanisms 



We discussed above the microscopic contributions to the relaxation rates arising from 
superconducting quasi-particles in the regime of Gaussian fluctuations. Now we proceed 
to study the macroscopic contributions of fluctuations due to VA excitations. The vortices 
cause contributions from both the quasi-particle excitations in their normal cores and from 
the fluctuating magnetic fields carried by them. These fluctuating fields interact directly 
with the nuclear magnetic moments and cause spin flip-transitions. The nuclear spins in the 
vortex cores and in the superconducting regions can be brought into thermal equilibrium by 
cross relaxation through simultaneous spin flips of neighboring nuclei, i.e., by spin diffusion. 
In type II superconductors, the cross relaxation of the nuclear spins in the normal core 
with the nuclei not in the core tends to be suppressed by the mismatch of the Zeeman 
energies in an inhomogeneous magnetic field.0 However, in high-T c superconductors in the 
absence of an external field cross relaxation is apparently not suppressed by this effect. The 
experimentally observed line width, of the order of T 2 _1 , is about 10 5 s _1 . On the other hand, 
the magnetic field variations due to the vortices and the corresponding variation of Zeeman 
energies of the 63 Cu nuclei lead to an inhomogeneous line width which is much smaller, of 
the order of 10 3 s _1 .il Hence the system is homogeneous; the 63 Cu nuclei can be considered 
as a single system subject to different relaxation mechanisms. Furthermore, the motion of 
the vortex cores, which visit many nuclei, also homogenizes the spin system. 

We now focus on the relaxation of nuclear spins caused by the fluctuating magnetic 
fields of the VA pairs. Since a static magnetic field is not easily taken into account in BKT 
theory, effects of BKT vortex fluctuations are best studied in the relaxation rates measured 
in nuclear quadrupole resonance (NQR) experiments on 63 Cu nuclei which have spin 1 = 3/2 
and nuclear quadrupole moment Q = —0.157 • 10~ 24 cm 2 . The local electric field gradient 
is oriented in the z direction and gives a quadrupolar splitting that is much larger than the 
Zeeman splitting due to the magnetic field of the vortices. The relaxation rates are governed 
by the time dependent correlation functions of the magnetic field, which are determined by 
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the VA fluctuations. The evaluation of these correlation functions for the diffusing and 
recombining vortices is the main task. We then proceed to discuss the NQR relaxation rates 
in terms of the magnetic-field fluctuations. 



C. Magnetic-Field Correlation Functions 

The Redfield theory gives the nuclear spin relaxation rates in terms of correlation func- 
tions of the fluctuating magnetic field.ilS We are interested in the correlation functions for 
the time-dependent local magnetic field in a given layer, 



k a p(t) = h nt a(r, t)h nt p(r, 0) = h OtCe (0, t)h 0t p(0, 0), (36) 

where h n a (r,t) is the a component of the magnetic field at point r and time t in layer n. 
This field originates from the VA pairs in all the layers and is given by 

N 

ho(0,t)=EE[ H -n(- r -n,,+ W) 

n u=l 

-H_ n (-r_ njl/ _(t))], (37) 

where r n v+ (t) and r n v -(t) are the positions of the vortex and antivortex of the i^-th pair in 
layer n at time t, and H n (r) is the magnetic field in layer n at r of a single vortex centered 
at the origin. The components of this field parallel and perpendicular to the layers are0'0 

H Hl( r ) = 4ZT^sign(n) 

ab 
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^ 4^r» + n V —) ■ <39) 

Here 4> = hc/2e is the flux quantum, X ab is the penetration depth inside the layer, and s is 
the inter-layer separation. Whereas the z component is of short range in both the in-plane 
and the z direction, the in-plane component is of short range only in the z direction and 
falls off with 1/r in the plane. For the convenient evaluation of the relaxation rates we will 
later use the Fourier transforms 
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for n ^ 0, H ||(k) = 0, and 

H " M = £k «p(-M ■ («) 



u \/* 2 + ^ 

We take the correlations between the fields of the vortex and the antivortex of the same 
pair into account but neglect inter-pair correlations. This is a good approximation, since 
the typical pair size is small compared with the average distance between pairs below the 
transition and even in a significant temperature range above it.EH In the following we are 
only interested in the diagonal components of k a p. They can be written as 

kaa(t) = / d 2 r' + d 2 r'__d 2 r + d 2 r- H^ a (-r' + ) 

? n J 

x [#_ n>Q (-r + ) - F_ n , a (-r_) 

xP(r / + ,r / _;r + ,r_;*)/(r + -r_). (42) 

Here, we assume the presence of N vortices and N antivortices in each layer so that N/ F = n 
is the vortex density calculated as a function of temperature in Refs. [12] and |38|. The function 
/(r) gives the normalized size distribution of the pairs. The diffusion function P describes 
the motion of the pairs within a layer. We now discuss the functions P and /. 

We assume diffusive motion of vortices but take the intra-pair interaction into account. 
The diffusion function P is defined as follows: P (V, , r'_; r + , r_; t)d 2 r' + d 2 r'_ is the probability 
of finding the vortex of a given pair in the area d 2 r' + about r' + and the antivortex of the 
same pair in d 2 r'_ about r'_ at time t, provided the vortex was at r + and the antivortex at r_ 
at t — 0. The diffusion function P is discussed in a recent paper by TimmE! in the context 
of flux noise. Here we summarize the results. Assuming for the moment that vortices and 
antivortices are unbound, P is simply a product of free diffusion functions for r + and r_ 
with a diffusion constant D. It can be rewritten as 
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1 / |R'-R| 2 \ 

"2vr J Dt eXP l v 2Dt ) 
1 / |r' -r| 2 \ 

X 8riS eXP ("^rj' <43) 

which is the product of free diffusion functions for the center of mass R and the separation 

vector r of a pair, showing that the center of mass and the separation vector diffuse with 

D cm — D/2 and D m \ = 2D, respectively. The assumption of unbound pairs would be justified 

for a small density of essentially free vortices, a situation that does not arise in practice. 

We take into account the interaction between vortex and antivortex, which is logarithmic in 

distance r, V(r) ~ q 2 ln(r/r ), where q is the charge of the vortex in the Coulomb gas model 

and r can be chosen as the size of the vortex core, and solve the Fokker-Planck diffusion 

equation. The result is 

P ( R' + -, R' - -; R+ -, R - -; t ) 

V 2' 2' 2' 2' ) 



IR' - R 



2 N 



=2^ ex 'v^^r' d(r '' r;() ' < 44) 

where the first term accounts for the motion of the center of mass and P re i gives the relative 
motion in polar coordinates r, 0, 

1 fry ( r' 2 + r 2 \ 
Prel ~ 4irD Tel t [r) eXP { AD Iel t ) 

"JL^'M^Y (45) 

Here, 7 = (1 — q 2 /ksT)/2 and I^{x) is a modified Bessel function. Note that P incorporates 
the effect of pair recombination: Pairs with zero separation are taken out of the process. 
Thus P starts out normalized to unity at t — but then drops to zero on the time scale 
of the life time r r of a VA pair, determined by D m \ = 2D. Newly created VA pairs do not 
enter in k a p, since there positions are assumed to be uncorrelated to those of existing pairs. 
If the diffusion of a vortex were limited only by the Bardeen-Stephen friction mechanism@ 
(no pinning) with a large diffusion constant, r r would be much smaller than the time scale 
of NQR spectroscopy, of order 10~ 7 to 10~ 8 s. We will come back to this problem below. 
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Besides the diffusion function, the correlation function in Eq. fl42|) contains the distri- 
bution function of pair sizes, f(r). Taking into account that the magnetic field of a vortex 
changes on the scale of the penetration depth X a b, and that for this reason the fields of a 
vortex and an antivortex almost cancel for separations r <C A a b, we approximate the pair 
distribution function with an analytical expression that becomes exact for large pairs and 
does not introduce irrelevant complications for small r, 

'<'> « T^T^. W> 

which should be normalized to unity. The form of the exponent ( is given in Ref. We 
only note that ( vanishes for T > T c and is positive and, to leading order, proportional to 
(T c — T) 1 / 2 below T c , where we now denote the BKT transition temperature by T c . 

D. NQR Relaxation Rates 

We can now evaluate the NQR relaxation rates in terms of the correlation functions, 
Eq. (fE|). The NQR rates are of interest here because the BKT vortex fluctuations show up 
more clearly in the absence of Abrikosov vortices due to an external magnetic field. In a 
vanishing field, however, the occupation of the quadrupolar energy levels cannot be described 
in terms of a Boltzmann distribution with a spin temperature and, therefore, we proceed by 
calculating the NQR relaxation rates using the Bloch-Wangsness-Redfield theory.0 We can 
apply this theory, since the interaction between the nuclear spins and the magnetic field of 
the vortices is a small perturbation compared with the quadrupolar splitting. The energy 
levels of in-plane 63 Cu nuclei due to quadrupole splitting are shown in Fig. 8. Since the field 
gradient is oriented along the z direction and I = 3/2, the relaxation rates are given by 



Tf 1 = 9k xx (0) - 7k xx (oo) , (47) 



T, 1 = 



1 3 

-k xx {u) + -k xx (0) + k zz (0) 



(48) 



where 7„ = 7.1 G : s 1 is the gyromagnetic ratio of the 63 Cu nuclei and u = 1^+3/2—^+1/2 1 = 
|a;_3/2 — 6c?_i/2 1 - The correlation function given by the temporal Fourier transform of Eq. (^2| 
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has the form 

F -D re l J n ^ 

00 /"OO 

x J2 Jm(kr/2) / dr'r' l+ ~<J m (kr'/2) 

m=l, odd 

xRe l^-, (^V4 + ^/Aeir<) tf/^j (y/kyA + iu/D^r^ , (49) 

where J, i^, and J are Bessel functions, a = x,y,z, r< = min(r, r'), and r > = max(r, r'). 
This equation shows how the spectral density of the magnetic field fluctuations determines 
the relaxation rates, Eqs. (f47D and (|48|) . With the Fourier transforms of the vortex magnetic 
field substituted from Eqs. (fK]) and (|4T|) , the final form of the correlation functions is 



*_(„) = *1 f * * /(r) 

teA al' ^ fc exp 2 SV ^T1^) - 1 Jo 

x £ J m (kr/2) / dr'r' l ^J m {kr'/2) 

1 jj Jo 
m—l, odd 



xRe + io;/Aeir< j ^y^i (y*74 + ^/Aeir> j (50) 

and k zz (uj) follows with a similar form. 

The effect of motional narrowing is taken into account: The nuclei are fixed and expe- 
rience the field fluctuations of the vortices as they pass by. An important parameter for 
diffusive relaxation phenomena is the characteristic frequency uo c corresponding to the in- 
verse jump time for the diffusion processEl Here, u c ~ D rc \/A\ 2 ah , the reason being that the 
dependence on to is determined by k 2 / '4 + iu j D re \ and the characteristic value of k is 1/A &. 
By inspection of k aa (u) it is seen that the correlation function becomes really frequency 
dependent only for frequencies uj > u c . The Fourier transform k aa (t), Eq. fl36|), begins to 
decrease for t < 1/uj c . However, k aa (t) does not have the simple exponential decay form 
characteristic for the case of nuclei diffusing in a static, random magnetic field. 

For the special case of unbound pairs, see Eq. fl43|) , the result for T{~ takes the simple 
form0 

x _ 3 7 > P ^ D k 
1 1 - 72 . .2\2 / 



4tt cu 2 A 2 ff Jo 1 + ZW/u; 2 
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(51) 



where rip is the temperature-dependent pair density, A e ff = 2X 2 ab / s is the effective penetration 
depthjlfl and ^ is the in-plane coherence length. In the general case there is no such simple 
form. We now proceed to discuss the results in terms of the diffusion constant D. 



The most important parameter in Eq. (|50|) is the diffusion constant D re \ = 2D that 
governs both the free motion of independent vortices and the motions of a vortex and an 
antivortex towards each other, leading eventually to their recombination. In the absence of 
pinning, D is given by the Bardeen-Stephen formulaB 



where p n is the normal-state resistivity and d is the layer thickness. For Bi-2212, D is of the 
order of 1 cm 2 /s at low temperatures.il This value is so large that vortices and antivortices 
would recombine so fast that there is no time left for a large number of nuclear relaxation 
processes to occur. Futhermore, also when ignoring recombination processes, the free motion 
of vortices and antivortices is so fast that the rapidly changing magnetic fields at the nuclei 
lead to slow relaxation, similar to the case of motional narrowing. 

Slow diffusion rates and correspondingly long VA recombination times are crucial for 
obtaining measurable contributions to Tf 1 and T% . In real high-T c superconductors there 
are always inhomogeneities, i.e., doping defects or intrinsic crystalline defects (twin bound- 
aries, grain boundaries), which can pin vortices by their interaction with the normal vortex 
cores. A measure of the strength of the pinning potential is the energy E p for the thermally 
activated motion of a vortex. Pinning leads to an Arrhenius-type temperature dependence, 
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(52) 




(53) 
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In terms of D the mean time between flights, or the time of stay, is given by t p = £l s /4D 
, where £ e g is a measure of the flight distance.El Experimentally, D has been determined 
for different thicknesses of layered CuC>2 systems. For a one-unit-cell thick film of Bi- 
2212 Rogers et a/.S measure the low-frequency flux noise near T c . By analyzing the 
frequency dependence of in terms of diffusion noise, the authors determine a characteristic 
frequency uj c = D(T)/2(r) 2 , above which S<p oc tu~ 3 ^ 2 ; (r) is the average pair size. From 
these experiments the authors determine D(T), Eq. (^), and find a temperature-dependent 
activation energy for a single CuC>2 layer, E P (T) ~ E (l — T/T c0 ), with E^jks ~ 800 K. For 
more than one layer thick epitaxial blocks of DyBa2Cu307_ x , the activation energy E p is 
proportional to the number of layers, N, and begins to saturate at N = 3 to iV = 4.0 This 
is taken as direct evidence that the pancakes in two adjacent layers are coupled due to the 
Josephson effect and move as entities. With N up to 10 5 , large activation energies E p ~ 8 eV 
(consisting of the nucleation energy plus the pinning energy) are observed in epitaxial films 
of Bi-2212, where a BKT transition is experimentally found in the temperature dependence 
of the penetration depth.0 In such crystals, the BKT transition may be driven by thermally 
created pairs of VA lines through several layers. 

To sum up, the diffusion constant D for a pancake vortex in a single clean CuC>2 layer 
is large, D ~ 1 cm 2 /s so that the VA life time is short. However, in the real multi-layer 
systems where the BKT transition is observed,i3I! D can be small and thereby enhance the 
life time of a VA pair so that the relaxation of nuclear spins can accompany the translational 
diffusion of these pairs. The actual value of D can vary over many orders of magnitude; for 
this reason we treat D as an open parameter in the following discussion. 

Let us first comment on Tj - 1 ignoring the interaction and recombination of VA pairs, 
cf. Eq. (pTp. We assume a multilayer structure of alternating single layers of YBa 2 Cu30 7 _ a . 
and PrBa 2 Cu 3 7 _ :r ./; the distance between the superconducting Cu0 2 layers is s = 24 A. 
Using A a 6 ~ 1400 A, £ a 6 « 12 A, and taking the experimental value D = 2 ■ 10~ 4 cm 2 /s, 
measured by Fiory et aS, for a thin film of YBa2Cu307_ x near T c , we get from Eq. flSTp a 
rate of Xf 1 « 572 s _1 , which is comparable with the experimental values. 
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Next, let us take into account the interaction and recombination of vortices and an- 
tivortices. The time it takes for a vortex to move to its nearest antivortex depends 
on the VA separation and on D for the single-vortex motion. Since D is not known, 
we evaluate the relaxation rates Tf and T 2 _1 as functions of the diffusion constant of 
the separation vector, D Te \ = 2D. The curves shown in Fig. 9 are obtained by nu- 
merical integrations from Eqs. (|47|)-([49l) using parameter values that apply to Bi-2212. 
The magnetic penetration depth and the Ginzburg-Landau coherence length are given by 
A a6 (T)/A a6 (0) = U(T)/U(0) = \Jt c q/ (T c0 — T), where A a6 (0) « 2000 A, £ o6 (0) « 21.5 A, 
and T c o ~ 86.8 K is the mean-field transition temperature. The distance between lay- 
ers is s ~ 15.5 A. The parameter 7 = (1 — g 2 //csT)/2 depends on the VA interaction, 
q 2 = q l(T c - T), where k B T c /ql « 0.2053 and for T c we take the value 84.7 K.HI Figure 9 
shows that the prefactor 1/-D re i in Eq. ( f49|) dominates the dependence of both relaxation 
rates on -D re i; both rates fall off approximately as 1/D TC \. The reason is that a large value 
of D re \ leads to fast recombination.0 Then, the time-dependent correlation functions k aa (t) 
are narrow and their Fourier transforms k aa (uj), Eq. (^), are broad and quite small even 
at uj w 0. Hence, small diffusion constants are necessary to observe contributions from VA 
fluctuations to the relaxation rates. 

The diffusion constant has an exponential temperature dependence if there are pinning 
effects. The corresponding temperature dependences of T{ 1 and T^ x are evaluated with 
D from Eq. fl53"|) , using for E p (T) the experimental results of Rogers et a/.H The bare 
diffusion constant Do is given by Eq. (0), it depends on T through £ a b and p n , the normal 
resistance of a single layer. The experimental value of D(T C ) is large, not much smaller 
than 1 cm 2 /s. Assuming the experimental values of D(T) obtained by Rogers et al, the 
calculated relaxation rates in the vicinity of T c = 84.7 K are shown in Figs. 10 and 11. 
It is found that the temperature dependences of both rates are qualitatively similar to the 
quasi-particle contributions above and below T c . There is a sharp drop of the rates below T c . 
The vortex contribution to T 2 _1 is smaller than the contribution to T{~ . Since the observed 
spin-spin relaxation rates are larger than the spin-lattice relaxation rates, the latter are more 
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suitable for an experimental test. 

The crucial point in our example is the small absolute values of the calculated relaxation 
rates. The numbers obtained by assuming a large value of D are several orders of magnitude 
smaller than the experimental values. These are of the order of 10 2 to 10 3 s _1 above T c , 
where the relaxation is predominantly caused by quasi-particle excitations. Furthermore, we 
have mentioned above that quasi-particle excitations in the cores will also contribute to the 
relaxation effects from vortices. In a simple picture the vortex cores can be considered as 
normal-conducting regions, where normal-state-like excitations contribute to T{~ according 
to the Korringa law (TiT) -1 = Kq; here Kq is the Korringa constant in the normal state. 
In the superconducting state the vortex-core contribution to Tf 1 is approximately given by 

— !— = 2n P (T) < 6 (T) K , (54) 

J- l,core-' 

where 2np is the areal density of vortices and 7r£ 2 6 is the area of one core. The calculated 
temperature dependence of T{~ core for Bi-2212 is shown in Fig. 12. A very similar curve 
would result for T 2 ~ c 1 ore . The core contribution is much larger than the contribution from the 
VA magnetic-field fluctuations, Fig. 10. 

There are, however, several effects that can make an experimental observation of the effect 
of the magnetic-field fluctuations feasible. We have already discussed possible mechanisms 
by which the diffusion constant D can become smaller, namely pinning and coherent motion 
of stacks of vortices. Even if D is sufficiently small, however, the VA contribution to T-f 1 
and T 2 _1 is at best comparable with the quasi-particle contributions. In order to evade 
the in-plane quasi-particle effects, inter-plane ions can be used for NQR spectroscopy. The 
relaxation of inter-plane nuclear moments can be affected by the magnetic-field fluctuations 
originating from vortices and antivortices in the CuC>2 planes.S A possible candidate is 201 Hg 
in HgBa 2 CaCu 2 6+ 5 with the 201 Hg ion with spin / = 3/2 and abundance 13.2% residing 
between two CuC>2 planes. So far NMR spectroscopy has been carried out only on the 199 Hg 
nucleus with / = 1/2; an NMR spin-lattice relaxation time of T\ = 32 ms is found at room 
temperature.0 Although this value will increase with decreasing temperature if Korringa's 
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relation applies, the relaxation rate Ti 1 at T c could still be an order of magnitude larger 
than the contribution from magnetic-field fluctuations. Other inter-planar candidates may 
be 87 Sr with 1 = 9/2 and abundance 7% and 209 Bi with 1 = 9/2 and abundance 100%, both 
in Bi-2212. The isotope 43 Ca (J = 7/2) is probably not a good candidate, since its abundance 
is rather small. To extend our approach to inter-plane ions, the sums over squared Fourier 
transforms of the vortex magnetic field in Eq. ( f49|) have to be taken over the appropriate 
fractional values for n. There are other inter-plane ions with I > 1/2 that are candidates 
for NQR relaxation measurements in order to observe the magnetic fields caused by the 
diffusional motion of BKT vortices. A pertinent question concerns the relative magnitudes 
of the magnetic fields at the inter-layer sites that originate from VA fluctuations and from 
the fluctuating Cu moments in the Cu02 planes, respectively. 

IV. SUMMARY, CONCLUSIONS, AND COMMENTS ON EXPERIMENTS 

In this paper we studied the effect of 2D superconducting fluctuations on the nuclear 
spin relaxation of layered superconductors above the superconducting transition temper- 
ature Tc and outside the region where Josephson coupling between the layers leads to a 
narrow 3D critical region surrounding Tc- Hence we model the superconductor as a set of 
2D layers without interlayer coupling mechanisms that transfer Cooper pairs between ad- 
jacent layers. The fluctuations in each layer consist of the usual Gaussian fluctuations and 
the topological Berezinskii-Kosterlitz-ThoulessS (BKT) vortex-antivortex fluctuations. The 
Gaussian fluctuations are considered as the long wavelength amplitude fluctuations of the 
superconducting order parameter and the BKT fluctuations are taken into account as the 
shorter wavelength phase fluctuations of each layer. These two fluctuation effects dominate 
the superconducting fluctuation behavior in the 2D fluctuation regime (cf. Fig. 1) of layered 
systems. 

We first evaluated in Sec. II, for s-wave superconductors, the effect of amplitude fluctua- 
tions on the hyperfine relaxation rate T{ 1 for pure systems, with and without pair breaking 
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effects due to inelastic scattering. The amplitude fiuctuations affect the Pauli spin suscep- 
tibility x °f the itinerant charge carriers and thereby, via the hyperfine coupling, the decay 
of the nuclear polarization due to the spin-lattice relaxation. Our results for the relative 
effect of fiuctuations on Tjf in zero magnetic field depends strongly on the pair breaking 
parameter ate = h/^r^ksTc as is seen from Figs. 6 and 7. The results for clean systems 
without pair breaking are given by Eqs. - fl26f) and with pair breaking by Eqs. (|33|) 
- (^1). Assuming small pair breaking effects, «c ~ 0.01, and the temperature just outside 
the critical region above Tc, say T > l.OlTc, a fluctuation enhancement of Tf occurs that 
is of order 10 with a slow decrease as the temperature moves away from Tc, see Fig. 7. 

Let us comment on the feasibility of experimental verification of this enhancement effect 
and its temperature dependence. A difficulty can arise in observing this effect: The high 
frequency fields used in NMR experiments on metals, and especially superconductors, is 
shielded within the penetration depth of the radio-frequency field.0 For this reason most 
of the experiments are performed on powders with grain sizes smaller than the skin depth 
(in metals) or the penetration depth (in superconductors), or on thin films. In powders 
or thin films, however, the NMR lines are broadened by charge density fluctuations that 
emanate from crystal surfaces, an effect similar to the Friedel oscillations of the charge 
density surrounding an impurity in a simple metal. S This geometrical broadening effect will 
primarily affect the relaxation rate T 2 _1 for nuclear spin-spin coupling, i.e., the natural line 
width. The charge fluctuations can also affect T^ 1 because of the change of the electron 
density at Ep. Hence, the singular behavior of Tf near the transition, seen in Fig. 7, may 
be smeared out by surface effects.il Furthermore, the spread of c-axis orientations in powder 
samples can also smear the effect of fluctuations on T' 1 _1 .0 

At this point one must ask what the chances are of finding systems in which these prop- 
erties can be observed. Aside from possibily the electron-doped high-T^ superconductors 
mentioned in the introduction, there are other quasi-2D superconductors with s-wave pair- 
ing and BCS behavior. The relaxation time Ti of the 93 Nb nuclear spins in a single crystal 
of 2H-NbSe2 (with trigonal prismatic coordination of the Nb atoms) has been measured 
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by Wada.El The rate follows the Korringa relation in the normal state and increases expo- 
nentially in the superconducting state. Other dichalcogenides include the TaS 2 - :I Se x . layer 
cpmpounds which are BCS superconductors and can be intercalated with organic molecules 
to increase the interlayer separation to as much as 50 A.0 Finally there are the graphite in- 
tercalation compounds such as CgK which are BCS superconductors with rather low T^'sS 
One or the other of such quasi-2D systems may be a suitable candidate for observing the 
T\ fluctuation effects discussed in this paper. The experimentalist may prefer to measure 
T\ on a nucleus without a quadrupolar moment in order to avoid the line-broadening effect 
caused by electric field gradients near the metal surfaces of thin films or small particles. 

In Sec. Ill we studied the effect of Berezinski-Kosterlitz-Thoules vortex-antivortex fluc- 
tuations on the NQR spin-lattice and spin-spin relaxation rates, TJf and T% x , respectively. 
Here the symmetry of the superconducting order parameter will affect the field and current 
distributions of an individual vortex. However, for both s- and d-wave superconductors, the 
flux contained in a vortex is $o , the magnetic field not too close to the vortex core does 
not depend on the symmetry of the order parameter, and, therefore, this symmetry will not 
affect the magnetic field fluctuations caused by the translational motions of an ensemble of 
VA pairs. Hence, the results for the relaxation rates T{~ and T 2 _1 calculated in Sec. Ill 
can be applied to both s- and d-wave superconductors. One must keep in mind, however, 
that the relaxation processes of the VA pairs are not entirely the result of the translational 
diffusion processes. The relaxation processes due to quasi-particle excitations must also 
be taken into account and can be different for the two symmetries. We assume that clear 
experimental evidence exists for the BKT transition in quasi-two-dimensional systems, in 
particular for layered superconductors with weak interlayer coupling such as some high-Tc 
cuprates. We also assume that this coupling does not change the qualitative behavior of our 
results obtained for a stacked system of uncoupled layers. For uncoupled layers, the basic 
tenet of dynamical BKT behavior is that vortices and antivortices move diffusively with 
bonding and unbonding of pairs under the influence of random internal forces (pinning). 
We studied the effect of the fluctuating magnetic fields accompanying the motion of vortices 
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and antivortices on the longitudinal and transverse NQR relaxation rates T^f 1 and T 2 _1 in 
the vicinity of T C - Our procedure uses a Coulomb gas description of thermally created VA 
pairs. The NQR relaxation rates were calculated in terms of the time-dependent correlation 
functions of the fluctuating magnetic field of the vortices and antivortices, Eq. d^9|). The 
interaction and recombination of VA pairs were taken into account and they can drastically 
reduce the correlation time. From the results obtained for Tf" 1 and T 2 _1 , Eqs. (fI?D and ([181). 
it is seen that the vortex diffusion constant D(T) plays a crucial role. The reason is that 
D(T) determines both the translational motion of the vortices and also the recombination 
time of the VA pairs. From Figs. 9-11 it is seen that D must be sufficiently small in order 
to get contributions to the relaxation rates that are comparable to those caused by quasi- 
particle excitations. We discuss possible candidates for the experimental NMR spectroscopy 
of dynamical BKT effects in Sec. III. As for the interlayer quasi-particle contribution to 
the relaxation rates, one way to eliminate these contributions is to observe NQR on in- 
terlayer nuclei. Possible candidates are 87 SR and 209 Bi in Bi-2212 and 201 Hg in mercury 
cuprates. Also for NQR experiments on small particles or thin films, the surface effects 
can be important because the charge density oscillations induced by the surface will set up 
electric field gradients that can interact with the quadrupolar moments of the nuclei. Hence, 
an unambiguous observation of 2D fluctuation effects in NMR or NQR requires thoughtful 
experimental set-ups. 
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FIGURE CAPTIONS 

1. Temperature regions. The mean field transition temperature of a 2D layer is given by 
Tco- The regime of long wavelength Gaussian fluctuations with 2D character lies outside the 
shaded region of strong 3D fluctuations. Vortex-antivortex fluctuations exist above Tbkt in 
a single layer. In a 3D system of weakly coupled layers, coming from high temperatures, 
the 2D fluctuations are cut off at the cross-over temperature T cr , before the true transition 
temperature Tc is reached. 

2. Transverse susceptibility diagrams through first order in the pair fluctuation propa- 
gator, V. 

3. The three fluctuation contributions to Imx(w) = Z) q Imx(q,u;) for a clean system 



without pair breaking measured in units of 



2 

a \ hui Nq 



to J E F 2 



4. The fluctuation propagator Dasa function of momentum x = k/kp for several values 
of the pair breaking parameter ac- 



(ft) 



2 hw Np 
E F 2 



as a func- 



5. The fluctuation contribution Imx^ r (u;) measured in units of 
tion of the pair breaking parameter ac for several values of the reduced external frequency 
f2 = huj/2Ep. For the experimental value = 5 x 10~ 8 the effect of even very little pair 
breaking is drastic. 

6. The fluctuation contributions to [l/(TiT , )]pj j , Imxr(^) = im XFF,r + ^- m XFF,r + 



ImxQQ T i measured in units of J^-^ 



as a function of the pair breaking parameter 
a c for T = 1.037b and U = hu/2E F = 5 x 1(T 8 . 

7. The fluctuation contribution to [l/(7i7 1 )]pL , Imxr(^), measured in units of 
(?o) 1^^2 as a f unc tion of Temperature for a range of the pair breaking parameter ac = 

0.01, 0.02, 0.04, 0.06, 0.08, and 0.10 (from top to bottom) and Vt = hu/2E F = 5 x 10~ 8 . 
The details for small and large pair breaking are shown in (b) and (c). 

8. Energy levels of 63 Cu nuclei in the Cu0 2 layer due to quadrupolar splitting. The 
allowed transitions are indicated by arrows. 

9. The contributions of vortex-antivortex magnetic-field fluctuations to the longitudinal 
and transverse relaxation rates T-f 1 and T 2 _1 as a function of the diffusion constant D re \ = 2D. 
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The dotted line gives the contribution from k xx (u = 0) to T-f , corresponding to the first 
term on the right hand side of Eq. (f47|). The parameters used in Eqs. ([47j)-(|49|) are given in 
the text. 

10. Magnetic-field contribution to the longitudinal relaxation rate T-j -1 near T c for Bi- 
2212. 

11. Magnetic-field contribution to the transverse relaxation rate T 2 _1 near T c for Bi-2212. 

12. Contribution to the longitudinal relaxation rate Tf" from vortex core relaxation. 
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